#include "laminarFlow.cpp"

using namespace std;
using namespace cfd;

int main() {
  LaminarFlow<double> lF;

//  for (int i=0; i<10; i++) {
//    lF.oneStep();
  lF.plot();
//  }
  
/*    vtkDataSetReader *reader = vtkDataSetReader::New();
    reader->SetFileName("data.vtk");
    reader->Update();*/
 
  return 0;
}
